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Abstract 

Likely candidates for the global potential energy minima of (H20)n 
clusters with n < 21 on the (OOOl)-surface of graphite are found us- 
ing basin-hopping global optimization. The potential energy surfaces 
are constructed using the TIP4P intermolecular potentials for the wa- 
ter molecules (the TIP3P is also explored as a secondary choice), a 
Lennard-Jones water-graphite potential, and a water-graphite polar- 
ization potential that is built from classical electrostatic image meth- 
ods and takes into account both the perpendicular and parallel electric 
polarizations of graphite. This potential energy surface produces a 
rather hydrophobic water-graphite interaction. As a consequence, the 
water component of the lowest graphite-(H20)„ minima is quite closely 
related to low-lying minima of the corresponding TIP4P (H20)n clus- 
ters. In about half of the cases the geometrical substructure of the wa- 
ter molecules in the graphite-(H20)n global minimum coincides with 
that of the corresponding free water cluster. Exceptions occur when 
the interaction with graphite induces a change in geometry. A com- 
parison of our results with available theoretical and experimental data 
is performed. 

*Correponding author. E-mail address: jmgomez@ull.es 
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1 Introduction 



The interaction of carbonaceous materials such as fullerenes, carbon nan- 
otubes and graphite with atoms and molecules share many properties. In 
particular, in this work we will be concerned with the interaction between 
water and graphite. A deep understanding of the features and properties of 
this interaction is of particular interest in technological applications, such as 
those related with the use of water as a lubricant for graphite fT|!2], and, more 
indirectly, in the behavior of water at the nanometer scales when material 
related to graphite, such as carbon nanotubes, are present. Water-graphite 
interaction is also relevant in the design of corrosion-free combustion cham- 
bers and rocket nozzles, since water is a universal combustion product and 
graphite is an important surface material because of its chemical inertness 
under extreme conditions fS]. Other fields benefiting from this knowledge 
include the environmental sciences |4] and astrophysics [5], since graphite is 
a good candidate for the composition of nano particles and dust grains. 

Despite the natural abundance of water and graphite, relatively few exper- 
imental data are available for their interaction. Studies at low temperature 
(T = 85 K) and low coverage using temperature programed desorption and 
vibrational high resolution electron energy loss spectroscopy have shown that 
water is adsorbed non dissociatively on the graphite surface forming hydro- 
gen bonded aggregates with a two dimensional structure that changes into 
a three dimensional one upon warming [6j. The water arrangement for the 
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two-dimensional structure is unknown, as is also unknown the role played by 
small water clusters in the growth of these structures. 

Experimental information about the water-graphite binding energies and 
structural aspects is currently lacking even for the water-monomer adsorp- 
tion. To our knowledge, there are only the early water-graphite binding 
energy by Kieslev et al. (15.0 kJ/mol) [Tj and the more recent association 
energy reported by Kasemo et al. [6]. 

In the last few years some results from theoretical calculations have been 
made available. Some of these studies are concerned with macroscopic fea- 
tures of the water-graphite interaction. In this group we can include the work 
by Werder et al. [H], who fit an interaction potential form to experimental 
data for the contact angle of water nanodroplets on graphite surfaces. A 
similar scheme is used by Pertsin et al. to simulate lubricant properties from 
a water-graphite interaction that was previously fitted to ab-initio [9j and 
empirical data |10| . Finally Gatica et al. have used empirical water-graphite 
potentials to look for a wetting transition [TT] . 

Ab initio calculations have been recently reported. By using second-order 
Moller-Plesset perturbation theory. Feller et al. [T2] have provided the inter- 
action energy between a water molecule and acenes as large as C96H24; the 
value of 24 kJ/mol that was obtained for this energy seems to be unphysically 
high [Hj. This conclusion is confirmed by the recent theoretical calculations 
by Sudiarta and Geldart [H]. Using the same Moller-Plesset scheme for 
a water molecule on both hydrogen and fiuorine terminated acenes, these 
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authors demonstrated the important contribution of this boundary to the 
water-acene binding energy. After removing this effect and extrapolating to 
an infinite graphene they reported a value of 10.2 kJ/mol. Density func- 
tional theory (DFT) total energy calculations were performed by Cabrera 
Sanfelix et al. |T5] to study structural aspects of water layers on graphite. 
Absolute binding energies were not provided and no global energy minimiza- 
tion was done. These tasks are performed in later calculations using DFT 
tight-binding methods complemented with empirical van der Waals force 
corrections (DFTB-D) P, [16]. In these works, clusters with up to 6 water 
molecules on graphite are studied and the optimal structures and the bind- 
ing and association energies were provided. Besides, these results are com- 
pared with those obtained with integrated ONIOM (ab-intio B3LYP:DFTB- 
D-|-semiempirical PM3) methods [3]. In these studies graphite is represented 
by up to three-layer acenes. In all these DFT studies the acene boundary 
effects are not removed. 

A full empirical approach has been followed by Karapetian et al. [I3] 
by using the Dang-Chang model for the water-water interaction and a po- 
larizable potential model for the water-graphene interaction that includes 
dispersion-repulsion contributions by means of Lennard- Jones pairwise in- 
teractions. In this potential the polarization term is built by associating an 
isotropic polarizable center with each carbon atom. The interaction between 
these centers, when polarized, is neglected. Again, only clusters with up to 
6 water molecules are considered. 



4 



A similar empirical approach shall be followed in our present work. Be- 
cause of its ability to reproduce the structure of water clusters, we will use 
the TIP4P model for the water- water interaction [17]. In order to analyze 
dependence of our results on this choice, we shall also consider the TIP3P 
water- water interaction model ^7\. The water-graphite interaction analytic 
model shall include a dispersion-repulsion term, built from the sum of infinite 
Lennard- Jones pairwise interactions using the Steele method, and a polar- 
ization contribution which shall be built using electrostatic image methods 
that take into account the anisotropic response of graphite. We shall provide 
likely candidates for the global potential energy minima of graphite-(H20)„ 
clusters with n < 21. We shall employ basin-hopping global optimization to 
identify these global minima. From the structure and energetics of these min- 
ima we shall elucidate about the hydrophobic nature of the water-graphite 
interaction at the lowest temperatures, and the dependence of the results on 
the potential model. 

This paper is organized as follows. In Section [2] we discuss our expression 
for the potential energy surface as a sum of Coulomb, dispersion-repulsion, 
and polarization contributions. In Section [3] we present likely candidates for 
the cluster global potential energy minima together with their association 
and binding energies. Here we shall compare our values with the available 
data. Finally, Section [4] summarizes our conclusions. 
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2 The Potential Energy Function 

The closed-shell electronic structure of both graphite and water makes an em- 
pirical approach to the potential energy surface (PES) for the water-graphite 
and water-water interactions particularly attractive. We write the potential 
energy of a graphite-(H20)„ cluster as a sum of two contributions 

V=V^^ + Kvg, (1) 

where V^w is the sum of pairwise water-water interactions, and Kvg is the 
water-graphite term. For the water-water interaction we have chosen the 
TIP4P form as a primary choice, but we will also consider the TIP3P model. 
These models describe each water molecule as a rigid body with two positive 
charges on the hydrogen atoms and a balancing negative charge either close 
to the oxygen atom (TIP4P) or just at the oxygen atom (TIP3P), together 
with a dispersion-repulsion center on the oxygen atom. Hence, V^^ is a sum 
of pairwise additive Coulomb and Lennard- Jones terms. These models have 
been used in the study of homogeneous water clusters [IH [HI [20], water 
clusters containing metallic cations [211 [22], and water-Ceo clusters [23] . 
The water-graphite interaction is written as 

Kvg = Vdr + Vpol, (2) 

where Vdr is a sum of pairwise dispersion-repulsion terms between the oxy- 
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gen and the carbon atoms. Each of these terms is expressed as a Lennard- 
Jones potential, whose parameters were obtained using the standard Lorentz- 
Berthelot combination rules from the corresponding parameters for the oxygen- 
oxygen and carbon-carbon interactions in TIP4P and TIP3P water and Steele 
[24] graphene-graphene potentials, respectively. Specifically, we used the val- 
ues eco = 0.389 kJ/mol and aco = 3.28 A for the TIP4P, and eco = 0.385 
kJ/mol and aco = 3.28 A for the TIP3P, which are similar to those derived 
by Werder et al. |8] to fit the contact angle for a water droplet on a graphene 
surface. An analytic form for Vdr can be obtained using Steele summation 
method [24l [25] over the graphite periodic structure by writing the interac- 
tion, [/dr, of a dispersion center at the point {x,y,z) with a graphite layer 
located at the surface z = (the origin of the reference frame is chosen at 
the center of a carbon hexagon), as a Fourier series, i.e. 



Udrix,y,z) = Uo{z) + J2Ui{z)fi{x,y) 

l>0 



(3) 



We have checked that the contribution to this expansion from terms with 
/ > 1 is negligible. Up to / = 1, we have 
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and 



ri(x, w) = — 2 < cos — xH ^ + cos — \x ^ + cos ^ — >, 




modified Bessel function of rrr order, and fi{x,y) is the first corrugation 
function. The total dispersion-repulsion interaction is obtained as a sum of 
t/dr terms over each graphite layer. We have obtained well converged values 
by including the Uq contribution from the two upper layers and the first 
corrugation of the first layer. 

In Eq. (I2|), V^oi includes the energy associated with the polarization of 
graphite due to the electric field of all the water point charges. This many- 
body interaction, which will turn out to be smaller than Vdr, was evaluated 
using a continuum representation of graphite. We will evaluate two contri- 
butions to Vpo\, 



each one associated, respectively, with the response of graphite to the elec- 
tric field component parallel and perpendicular to the graphite surface. For 
the first one, V\\, we will assume that graphite behaves as a classical conduc- 
tor, which allows us to make use of the image charge method to obtain (in 
Gaussian units) 



V^poi = Ml + V^, 



(7) 




(8) 



8 



where Qi is each of the water electric point charges and r^^ is the distance 
between the charge and the image of the charge qj, i.e., r'^j = [{xi —Xj)^ + 

In order to evaluate V±, we will assume that a graphite layer (at z — 0) 
has, in reciprocal space, a surface polarizability density Q:±{kxky) such that 
when an electric field depending on the surface point and perpendicular 
to the layer, E±{x,y) — ^ J dkxdkye'^^'^''^'^''y''^^£±{kx,ky), is applied, an 
electric dipole density, I±{x,y) = j dkxdkye~^^^''^^^y'^^X±{kx,ky), is in- 
duced on that layer, with I±{kxky) — a±{kxky)S±{kx, ky). If we now ne- 
glect the dependence of a± on kx and ky (which is a valid approximation 
if E±{x,y) depends smoothly enough on the surface point), then we would 
have I±{x,y) = a±E±{x,y), with a± = a±{0,0) being the electric polariz- 
ability in a uniform electric field perpendicular to the layer. In this way we 
can calculate the dipole density induced in the graphite layer by an electric 
charge Qi at the point {xi, yi, Zi), and from this dipole density we can evaluate 
its electric field and the electric force between the polarized layer and that 
charge. One readily shows that the electric field due to the polarized surface 
in the half-space of the charge {z > 0) is equal to the electric field induced by 
an image dipole = —2Tra±qi at the point {xi,yi, —Zi) and direction parallel 
to the z axis. This result can be generalized additively to the case of several 
electric point charges, all of them located in the space region ^ > 0. From 
the corresponding image dipoles we can obtain their electric force on each 
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charge, and from here the interaction potential V±, namely 

^2TTa:iqf 1 ^27ra±qiqj{zi + Zj) 

^^ = -2^—^2 nl^ 33 ■ (9) 

i ^ ijtj 'ij 

The use of expression ([9]) requires the knowledge of a±. The value of this 
polarizability density may be estimated from e±, the relative electric per- 
mittivity of graphite for applied electric fields perpendicular to the (0001) 
surface, whose value is e± = 5.75. Following Hannay's alternative deriva- 
tion of the Clausius-Mossotti equation [26], we shall require the form of the 
diverging term in the expression of the electric field induced by a uniform 
surface dipole density of magnitude I± and direction perpendicular to the 
surface. The magnitude of this term is readily found to be —4:nl±6{z). By 
space averaging it, we obtain the local electric field on each graphite layer 
-^locai = E + 47rlj_/d, where d = 3.35 A is the layer-to-layer distance and 
E is the magnitude of the macroscopic electric field in the medium. Then 
from the relation I± = a±Eioca.i, and the expression P = ^'^^^^E relating the 
volume polarization density P = I±/d and E, we arrive at the desired result 

We obtain by this procedure the value a± = 0.220 A. 

Being consistent with our metallic assumption for V\\, in the evaluation 
of the polarization contribution to the water-graphite interaction potential 
we shall assume total screening of the electric field by the external graphite 
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surface. Therefore, only the most external graphite layer shall be considered 
in this evaluation. 

We have also checked the relevance of the McLachlan substrate mediated 
interaction [27] between the water molecules in the presence of the conducting 
graphite layer and found it to be negligible (~ 0.03% of the total interaction 
energy for the water dimer on graphite); therefore, we shall not include this 
term in our potential energy. 

As mentioned in the Introduction, Karapetian et al. [13] have proposed 
a different model potential for the polarization contribution to the water- 
graphite interaction energy. These authors locate an isotropic polarizable 
center at each carbon atom and calculate the polarization energy as a sum of 
the contributions from each center. This model does not take into account the 
collective properties of the delocalized vr electrons and neglects completely 
all screening effects among the induced dipole moments. We have estimated 
that these defects of this model lead to an overestimation of the polarization 
energy by a factor of two for the water monomer. 

3 Global Potential Energy Minima 

Likely candidates for the global potential energy minima of graphite-(H20)„ 
clusters with n < 21 were located using the basin-hopping scheme [28], which 
corresponds to the 'Monte Carlo plus energy minimization' approach of Li 
and Scheraga [29]. This method has been used successfully for both neutral 



11 



[28] and charged atomic and molecular clusters [23l [30l [3H [32l [33] , along with 
many other applications [34]. In the size range considered here the global 
optimization problem is relatively straightforward, but somewhat more costly 
than in the water-Ceo clusters [23] • The global minimum is generally found in 
fewer than 7 x 10^ basin-hopping steps, independent of the random starting 
geometry. In some cases, starting out from the (H20)„ global potential 
minimum, the corresponding global minimum for graphite-(H20)n is found 
even faster. 

For graphite-(H20)n clusters, association energies, Ai?a, are defined for 
the process 

graphite + nHaO = graphite- (H2O),,; -AE^. (11) 

We also define the water binding energy, Ai^b, as the difference between the 
association energies of graphite- (H20)n and (H20)„, i.e. 

graphite + (H20)„ = graphite-(H20)„; -AE^. (12) 

The clusters in these expressions are assumed to be in their global minimum. 
The structures and association energies employed here for the global minima 
of (H20)„ coincide precisely with those obtained by Wales and Hodges [T8] 
and Kabrede and Hentschke [T9] . 

For comparison with the available data, we plot in Fig. [T] the association 
and binding energies defined in (fTTll and (fT2ll as a function of the number of 
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Figure 1: Association, AEg,/n (a), and binding, AE\y (b), energies per water 
molecule, for the global minima of water-graphene clusters. Our results: 
circles for TIP4P and diamonds for TIP3P. Other calculations: DFTB-D 
p] (crosses); DFTB-D [3] (asterisks); empirical [l3] (up triangles); ONIOM 
[3] (down triangles); Moller-Plesset [O] (square, only for n = 1). 
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water molecules, n < 6, for water on graphene. Our energies for the TIP3P 
model are somewhat higher than those for the TIP4P. With the exception 
of the ONIOM data, our binding and association energies are systematically 
lower and higher, respectively, than the ab initio and other empirical poten- 
tial values. However, it is well known that DFT methods tend to overestimate 
the binding energies. Furthermore, the acene-boundary effects discussed in 
the Introduction have not been removed in the ab initio calculations, except 
in Sidarta and Geldart's binding energy for the water monomer. It is also 
relevant to remind here the overestimation of the polarization energy that 
takes place in the empirical potential model by Karapetian et al, as men- 
tioned at the end of Section [21 Taking all these facts into account, we can 
conclude that the binding energies provided by our interaction model are 
quite reliable. 

Associations energies for n > 2 are dominated by the water-water in- 
teraction. As far as the empirical potential data are concerned, differences 
in these association energies are a consequence of the different water-water 
model interaction used in each case. The lower ab initio values reported 
by Xu et al. can be attributed to the poor description of the water-water 
interaction in their DFTB-D scheme. 

The structures of the lowest minima obtained for graphene-(H20)„ with 
1 < n < 6 coincide with those provided by the other available calculations 
and with the structures that we will obtain later for the graphite-(H20)„ 
clusters. For n = 1, our calculation, as the one done by Karapetian et al, 
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provides an H2O molecule with an OH bond pointing towards the graphene 
and the oxygen atom just over the center of a carbon hexagon ring. On the 
other hand, the ab initio data provide a two-legged conformation with the 
two OH bonds pointing towards graphene in a symmetric way. This might 
be understood as the result of quantum zero-point-energy effects, since the 
difference in energy between the two geometries may be of the order of the 
zero point energy. 

The association {AEg,/n) and binding energies (AE'b) for the full graphite- 
(H20)„ clusters, which have been calculated as described in Section [2] with 
the TIP4P water-water interaction, are plotted in Fig. [2l We also in- 
clude in Fig. [2](a) the values of the polarization energy Vpo\ and water- 
graphite dispersion-repulsion energy Vdr, as defined in Section [2], for the clus- 
ter global minima. The term Vpoi oscillates with n around an average value 
of V^poi = 3.33kJ/mol; the two contributions to V^oi, V\\ and V±, are similar 
in magnitude with V|| somewhat larger than V±. The term Vdr fluctuates also 
around a slowly growing average as the number of water molecules close to 
the graphite surface increases. On average, each of these water molecules 
contributes about 7.26kJ/mol to Vdr- The water-graphite binding energies 
correspond quite closely to the sum of Vpoi and Vdr, while the association 
energies are dominated by the water-water interaction. The average value of 
the association energy per molecule in homogeneous TIP4P (H20)„ clusters 
with 6<?T,<21is~ 42kJ/mol [IHl IHj. For water cluster on graphite the 
corresponding value turns out to be 44.6kJ/mol, which is comparable with 
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Figure 2: (a) Polarization Vjpoi (dashed line) and dispersion-repulsion Vdr 
(full line) contributions to the potential energy of global minimum graphite- 
(H20)„ clusters, (b) The corresponding association energies per water 
molecule, /S.Ea_/n (full line), and binding energies, Ai^b (dashed line). Global 
minima in which the structure of the (H20)„ moiety differs from the global 
minimum of the corresponding TIP4P (H20)n cluster are marked according 
to the discussion in the text. 
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the experimental value of 43.4 ± 2.9kJ/mol [6]. Any of these values corre- 
sponds to the binding energy of a water molecule in a water cluster, and it is 
much larger than the energy for binding a water molecule onto the graphite 
surface. This energy balance would support an hydrophobic nature of the 
water-graphite interaction. 

For n = 1 we obtain a binding energy /S.E^ = AEy^ = 8.81kJ/mol. 
This value is somewhat larger than our TIP4P binding energy for a water 
molecule on graphene (Ai?b = 7.6kJ/mol), and for the corresponding Ceo- 
(H2O) cluster {AEy, = 6.31kJ/mol) [23j. Thus, the numbers provided above 
for the binding energies of these three compounds are, at least, physically 
consistent. In the present case, the contribution of the polarization energy to 
the graphite-H20 binding energy is 2.08kJ/mol; the corresponding value in 
the C6o-(H20) cluster was 2.32 kJ/mol [23]. This larger value is a consequence 
of the important small-size quantum effects that make the polarizability of 
the Ceo molecule significantly larger than that of a conducting sphere with 
the geometrical Ceo radius [35]. The polarization energy is responsible for 
orienting the H2O molecule with an OH bond pointing towards the graphite 
surface and the oxygen atom just over the center of a hexagonal carbon ring 

(Fig. ED. 

The angle between the water C2 symmetry axis and our ^-axis in the 
n = 1 global minimum is 39.8 degrees, practically identical to the corre- 
sponding value in the Ceo-(H20) cluster (40.4 degrees) [23] and close to the 
experimental value in benzene (37 degrees) [36]. A different water orientation 
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Figure 3: Two views of the global minimum obtained for graphite- (H2O). 
This figure, and those that follow, were prepared using the program XCrys- 



with the two OH bonds pointing towards the graphite surface is also a local 
minimum, but it has a slightly lower binding energy (Ai?b = 8.67kJ/mol). 
The energy difference between this minimum and the global minimum (~ 
0.14kJ/mol) is so small that zero point energy effects might as well favor 
the two-legged structure as the vibrationally averaged quantum global min- 
imum [H]. The equilibrium distance in the global minimum between the 
oxygen and the graphite surface is 3.12 A, which is very close to the ab initio 
value (3.04 A) [I^ and the corresponding values in water-Ceo (3.19 A) [23] 
and water-benzene (experimental, 3.33 A) |36| . 

The structures of the TIP4P lowest minima obtained for graphite- (H20)n 
are presented in Fig. [H 

Due to the hydrophobic nature of the water-graphite interaction, the wa- 
ter substructure is often very similar to that in the corresponding global 
minimum of TIP4P (H20)„ [HI [19]. In some cases the structures are ac- 
tually identical (aside from minor differences in angles and distances). The 
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exceptions are labeled in Fig. [2j for those indicated by an asterisk (n = 
6, 7, 11, 15, 17, 19, 21) the water substructure corresponds to a low-lying lo- 
cal minimum of TIP4P (H20)„, rather than the global minimum. The energy 
penalty for this choice is mainly compensated by a more favorable dispersion- 
repulsion contribution to the interaction energy with graphite, which arises 
from a larger water-graphite contact area. For example, in the graphite- 
(H20)6 global minimum the geometry of the water moiety corresponds to 
the "book" structure, which has also been identified as the lowest minimum 
in the corresponding water-Cgo cluster [23]. For the sizes in question, there 
also exists a higher energy local minimum in which the water substructure 
corresponds to the global minimum for TIP4P (H20)„. The difference be- 
tween the association energies of the global and local minimum structures of 
each of these graphite-(H20)„ clusters is 6.36, 4.11, 1.62, 1.66, 13.04, 16.39 
and 9.41kJ/mol for n = 6, 7, 11, 15, 17, 19, and 21, respectively. 

For the clusters labeled with a circle in Fig. [2] = 9, 12, 13) the water 
substructure corresponds to a perturbation of the TIP4P (H20)„ global min- 
imum. In other words, relaxing the water moiety in the absence of graphite 
does not lead to a nearby local minimum, in contrast to the cases above. 
These graphite-(H20)„ structures appear to be favored by polarization con- 
tribution in the case n = 12 (only the hydrogen-bond pattern is modified), 
and both dispersion and polarization contributions for the other two cases. 

Notice that the structure of the water moiety up to n = 6 are basically 
planar with an average oxygen-graphite distance zq = 3.22 A; the n = 7 
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cluster corresponds to a transition to two-layer water clusters (^o = 4.04 
A), and in the range 8 < n < 21, we have always two-layer water structures 
(zo = 4.53 A) in which the clusters with odd n have one more water molecule 
in the layer closer to the graphite surface. 

The complete two-layer water structures for even n are precisely the struc- 
tures of the global TIP4P free water clusters. Therefore, these structures 
interact with graphite in an optimal way and they keep their structure in 
the corresponding water-graphite clusters. On the other hand, for odd n, 
the free water global minima do not show optimal surfaces for its interaction 
with graphite, thus explaining why these clusters change their structure to 
minimize that interaction energy. The chosen new structures are sensibly 
determined by those of either the n — 1 or n + 1 clusters. Other water- water 
potential models do not produce this alternating behavior in the structure 
of the free water global minima and, therefore, we can expect also different 
behavior in the water- graphite global minima for n > 8, as some preliminary 
results with the TIP3P model seem to confirm. 

The alternating behavior found in the structures of the water-graphite 
global minima determines the behavior of the second energy differences. 
These account for the relative cluster stability and their values for associ- 
ation and binding energies, per water molecule, are plotted in Fig. [5l For 
n > 7, we observe an oscillation of period An = 2, which corresponds to the 
discussed even-odd alternating structures. The data for the association ener- 
gies indicate that clusters with even n are more stable than their neighbors. 



21 



4i 




2 4 6 8 10^12 14 16 18 20 



Figure 5: Second energy differences per water molecule for the association 
energies (full line) and binding energies (dashed line) of water-graphite clus- 
ters. 

This pattern is in complete correlation with the corresponding behavior of 
the TIP4P free water global minima. Fig. [5] shows that, in general, second 
differences for binding and association energies are anticorrelated. The in- 
creased stability seen in the binding energies of odd n clusters relative to 
their neighbors is due to their extra water molecule in close contact with 
the graphite surface. From these features we conclude that the water-water 
interaction certainly dominates the observed relative cluster stability. Par- 
ticularly stable clusters occur for n = 4, 8, 12, 16, 20. 

In the light of the preceding results, one could ask a couple of questions 
that are relevant to assert the hydrophobic nature of the water-graphite in- 
teraction: Are there for n > 7 water-graphite local minima structures, close 
in energy to the global minima, in which the water molecules grow into a sin- 
gle layer (wetting structures)? How much larger should the water-graphite 
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interaction be for the previous wetting structures to become the most sta- 
ble ones? Answering to the second question, we have found that we had to 
multiply the value of the Lennard- Jones parameter eco by almost a factor 
of two for a wetting structure to become the global minimum. This is con- 
sistent with the analysis performed by Werder et al. [8] for the monomer 
binding energy required to produce a wetting behavior. By relaxing the wet- 
ting structures found by this procedure to the closest local minimum of our 
original potential we have found that these wetting local minima lie ~ 2.2 
kJ/mol per water molecule above the global minimum. Although these val- 
ues are smaller than those found in C6o(H20)„ clusters [23], the hydrophobic 
nature of the water-graphite interaction is also a quite robust property that 
would require unphysical changes in our model potential to modify it. 

We have already shown by making use of the TIP3P potential that the 
structures of the global minima of the first six water-graphite clusters are 
going to depend weakly on the model chosen for the water-water interaction. 
However the dependence found in the structure of the water-graphite global 
minima on the structure of the corresponding free water clusters and the 
known dependence of the latter on the water-water interaction model for 
n > 6, would imply changes in the structure of these larger water-graphite 
clusters when a different water model is chosen. Preliminary results confirm 
this prediction, with changes that are more significant than those found for 
C6o(H20)„ clusters. 
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4 Conclusions 



Using a theoretically guided empirical potential energy surface and basin- 
hopping global optimization we have characterized the geometrical structures 
and energetics of likely candidates for the global potential energy minima of 
graphite- (H20)„ clusters up to n = 21. The structures of these minima for 
1 < n < 6 coincide with those provided by other available calculations. For 
n > 2, association energies are dominated by the water- water interaction 
while the main contribution to the binding energies comes from the disper- 
sion energy. Our potential energy surface provides a rather hydrophobic 
water-graphite interaction at the nanoscopic level. As a consequence of this 
property the water substructure in the lowest energy clusters often corre- 
sponds closely to a low-lying minimum of the appropriate (H20)„ cluster. 
In most cases the structure is simply a slightly relaxed version of the global 
minimum for (H20)„. However, the presence of graphite can induce changes 
in geometry of the water moiety. 

For n = 6, 7, 11, 15, 17, 19, 21 the water substructure is based on a local 
minimum of (H20)„, which is close in energy to the global minimum. The en- 
ergy penalty for this choice is mainly compensated by the dispersion-repulsion 
contribution to the interaction energy, because the change in structure gives 
rise to a larger water-graphite contact surface. For n — 12 the water sub- 
structure is based on a deformation of the free (H20)i2 with the same oxygen 
framework as the global minimum, but a different hydrogen-bonding pattern. 
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The different orientation of some of the OH bonds close to the graphite sur- 
face increases the polarization energy, which stabilizes the structure. Finally, 
for n = 9 and 13 the water substructure involves a more significant defor- 
mation of the global minimum, the energy penalty being compensated by 
both the polarization and dispersion terms. A clear alternating behavior in 
which the water moiety of the clusters with even n keep the structure of 
the corresponding free water clusters, while their odd n neighbors change it, 
has been observed; this behavior has been shown to be completely correlated 
with that of the free water clusters, namely (H20)„ clusters with even n > 6 
present faces that interact with graphite in an optimal way, while those with 
odd n do not show this feature. 

Our potential energy surface also supports, for n > 6, wetting local min- 
ima in which the water molecules grow into a single layer. The potential en- 
ergies of these structures lie at least 2.2 kJ/mol per water molecule above the 
global minima. The presence of the graphite surface is necessary to stabilize 
the monolayer which otherwise collapses in the cases we have considered. In 
this respect, we have also found that in order to make the wetting local min- 
ima to become the cluster global minima we should increase the magnitude of 
the graphite-water interaction to unphysically high values. This implies that 
the hydrophobic nature of the water-graphite interaction at very low tem- 
perature is a quite robust property that would require unphysical changes in 
our model potential to modify it. 

In order to study the dependence of our qualitative picture on the model 
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chosen for the water- water interaction (the TIP4P potential was our primary 
choice), we have repeated our calculations for the TIP3P potential for n < 
6. The global minimum structures found for these clusters coincide with 
those of the TIP4P model. However the dependence found in the structure 
of the water-graphite global minima on the structure of the corresponding 
free water clusters and the known dependence of the latter on the water- 
water interaction model for n > 6, would imply changes in the structure of 
these larger water-graphite clusters when a different water model is chosen. 
Preliminary results the TIP3P model potential confirm this prediction. 

The lowest energy structures obtained in the present work will be made 
available for download from the Cambridge Cluster Database [37] . 
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